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Abstract 

To understand the essence of the exciton Mott transition in three-dimensional electron-hole systems, the metal- 
insulator transition is studied for a two-band Hubbard model in infinite dimensions with interactions of electron- 
electron (hole- hole) repulsion U and electron- hole attraction —U'. By using the dynamical mean- field theory, the 
phase diagram in the U-U' plane is obtained (which is exact in infinite dimensions) assuming that electron-hole 
pairs do not condense. When both electron and hole bands are half-filled, two types of insulating states appear: the 
Mott-Hubbard insulator for U > U' and the biexciton-like insulator for U < U' . Even when away from half-filling, 
we find the phase transition between the exciton- or biexciton-like insulator and a metallic state. This transition 
can be assigned to the exciton Mott transition, whereas the Mott-Hubbard transition is absent. 

Key words: exciton Mott transition, clcctron-holc systems, two-band Hubbard model, infinite-dimensions, dynamical mean-field 
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1. Introduction 

Electron-hole (e-h) systems in photoexcited semi- 
conductors exhibit various remarkable properties de- 
pending on carrier density, temperature, etc., and have 
been investigated extensively both experimentally and 
theoretically [1]. In particular, the metal-insulator 
transitions have attracted interest for many years: the 
exciton Mott transition at high temperatures between 
an exciton/biexciton gas phase and an e-h plasma 
phase, and crossover at low temperatures between the 
Bose-Einstein condensation (BEC) of excitons at low 
density and the BCS-like condensation of e-h pairs at 
high density. However, the complicated tangle of the 
elements, two types of fermions. Coulomb interactions 
of both repulsion and attraction, screening efi'ects, e-h 
densities, temperatures, etc., makes the physics of this 
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system hard to be understood. Therefore, theoretical 
understanding of especially the exciton Mott transi- 
tion and the BEC-BCS crossover is still not sufficient. 

We expect that a study from a standpoint of the 
strong-correlation physics provides new interpretation 
about the above problems. As the first step of our work, 
we examine the exciton Mott transition in considera- 
tion of the minimum elements, i.e., a two-band Hub- 
bard model, by using the dynamical mean-field the- 
ory (DMFT) [2] recently developed through the study 
of strongly correlated electron systems. The DMFT 
requires only the locality of the self-energy, and can 
take full account of local correlations. This locality and 
the resulting DMFT become exact in the limit of infi- 
nite spatial dimensions and good approximation of the 
three-dimensional systems. 

In the present paper, we focus on the normal phase 
where the condensation of e-h pairs (i.e., exciton BEC 
and e-h BCS state) is not allowed. The calculation is 
performed at absolute zero temperature. 



Preprint submitted to EXCON '04 Proceedings 



30 June 2004 



2. Two-site dynamical mean-field theory 

We consider a electron-hole system described by the 
two-band Hubbard model given by 

<ij>a a— c,h ia,(y. 

+ Uj2d?r^d?,d-!d?i -U'J2 d^UU^ldl,, (1) 

where d'^l (d^l) denotes a creation operator of a con- 
duction band electron (a valence band hole) with spin 
= {T; i} at the i-th site. The quantities te (th) and 
Mc (jJb) are the transfer integral of the electrons (holes) 
between the neighboring sites and the chemical poten- 
tial measured from the center of the electron (hole) 
band, respectively. The on-site Coulomb interaction of 
the e-e (h-h) repulsion and that of the e-h attraction 
are expressed by U and —U', respectively. The local 
Green function for electrons or holes of the model (1) 
is defined by 
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where T,"{uj) is the self-energy of electrons (q = c) 
or holes (q = h), which is local, i.e., does not de- 
pend on the wave number, in the limit of infinite- 
dimensions. Wc use the semicircular density of states 
(DOS), pTAe) = ^/Ul-e'/{2^ti). 

Within the DMFT [2], the many-body problem of 
the lattice model (1) is mapped onto the problem of a 
single-site impurity embedded in an effective medium. 
The effective medium, which is dynamical and is repre- 
sented by the noninteracting impurity Green function 
Qoiuo) of an effective single- impurity Anderson model 
(SIAM), is determined from the self-consistency con- 
dition Go{ujy^ = He - taG"(w). The condition 
is also read as Gf,^p(Lo) = G"(uj). The interacting im- 
purity Green function of the effective SIAM, Gf^^iuj), 
should be calculated exactly such that effects of the in- 
teractions on the impurity site are fully included. Con- 
trary to the ordinary mean-field approaches, thus, in 
the DMFT scheme the local correlations and dynami- 
cal quantum fiuctuations are taken into full account. 

In order to extract a sketch of the phase diagram of 
the model (1), here we apply the two-site DMFT [3] 
which is simplified version of the DMFT. In the two- 
site DMFT, the effective medium Qa{uo) is represented 
approximately by only the fewest parameters, i.e., the 
effective SIAM consists of a single impurity and only 
a single bath sites. Since the essence of the DMFT 
concerning the local correlations still remains despite 
the bold approximation and simplification, it can be 
successful to provide the most correct critical point 



of the Mott-Hubbard transition and to describe the 
qualitative electronic properties [3]. 

For the model (1), the corresponding effective two- 
site SIAM is written as 

//imp = C« + Vc [dpC^ + h.C.) - HadPd^] 

<7.a 

+ UJ2dr^d'^dfd'i -U'J2 rfa^rfa^irf^, (3) 
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where the bath parameters Va and e" denote the hy- 
bridization between the impurity (d) and bath (c) sites, 
and the energy level of the bath site, respectively. The 
Green function of the effective medium (i.e., noninter- 
acting impurity Green function) becomes So ('^)~^ = 
oj + Ha- V^/{ui - e'^). In the two-site DMFT, the self- 
consistency condition is reduced to simpler equation 
by the following procedure [3]: the self-energy is ex- 
panded in the low-energy region, S™(a;) ~ -|- baOJ, 
and then the resulting local Green function (2) and im- 
purity Green function Gg„p(w)-i = gS{<^)~^ - S"(w) 
arc compared so as to coincide at high-energy region. 
Thereby, the self-consistency equation for Va is ob- 
tained as 



where 

Za = {I — ba)~^ = 
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is the quasiparticle weight which generally character- 
izes the Fermi liquid (metallic) states. On the other 
hand, the requirement that the particle densities of the 

original and impurity models must be equal, i.e., = 
n"^p, leads to the self-consistency condition for e", 



j dLjlmG°'(ij + iO+) = j dwImGLp(a;-FiO+).(6) 
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Consequently, the model (1) can be solved within the 
two-site DMFT by the following self-consistency cycle: 

(i) G^p(w) is directly calculated by the exact diago- 
nalization of the two-site SIAM (3) with e" and Va- 

(ii) By using E"(a;) = So ('^)"' - GLp(a;)-\ a new 
value of Va is determined from the condition (4) and 
Eq. (5). (in) By substituting S"(a;) for Eq. (2), a new 
value of e" is chosen so as to satisfy the condition (6) . 
This process (i)-(iii) is iterated until e™ and Va con- 
verge. 

The metal-insulator transition for the normal phase 
of the model (1) is discussed from behaviors of both 
the quasiparticle weight Za and the interacting DOS 
p"[ij) = -ImG"(cj-f i0+)/7r, with varying [/, U',th/te 
and also the particle density n (= n'^ = n^). 
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3. Phase diagram at half filling 

First, we concentrate on the special case where the 
both electron and hole bands are half-filled, i.e., n = 1. 
In this symmetric case, we can set //a = f7/2 — U' and 
e? = 0. 

For fh/fc = 1 (the mass of the hole is the same as 
that of the electron), the phase diagram on the plane 
of U' and U is shown in Fig. 1. There are three kinds 
of states: (I) metallic state, (II) Mott-Hubbard insu- 
lating state, and (III) biexciton-like insulating state. 
The second-order transitions between these states oc- 
cur on the solid curves. In the metallic state (I), Za 
has a finite value and there is finite DOS at the Fermi 
level (the quasiparticle coherent peak), i.e., p"(0) ^ 0. 
On the other hand, in the both insulating states (II) 
and (III), Za = and the coherent peak of the DOS 
disappears. However, the physical pictures of the insu- 
lating states (II) and (III) are quite different, as drawn 
schematically in Fig. 1: the state (II) is induced by the 
e-e (h-h) repulsion U on each electron and hole band, 
while the state (III) is realized by the e-li attraction 
U' on each site. The competition of these two states 
stabilizes the metallic state for U ^U' . 

We point out that the above results are equivalent 
to those obtained for the two-orbital repulsive Hub- 
bard model [4] because this model and our model (1) 
only at half-filling can be mapped onto each other by 
the attraction-repulsion transformation. Actually, the 
phase diagram of Fig. 1 is in good agreement with that 
of Ref. [4]. 

We also examine effects of the difference between the 
electron and hole masses. In Fig. 2, the phase diagram 
on the plane of U' and U is shown for t^/te = 0.5, where 
the hole is twice as heavy as the electron. A new state 
(IV) appears between states (I) and (II) , in which Z^ ^ 
but Z-^ = 0, i.e., the electron (hole) band is metallic 
(insulating). In other words, the Mott-Hubbaxd tran- 
sition of holes does not coincide with that of electrons 




Fig. 1. Phase diagram in the U'-U plane at half-filling (n = 1) 
for th/te = 1. 
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Fig. 2. Phase diagram in the U'-U plane at half-filling (n = 1) 
for th/te = 0.5. 

when tc 7^ th- This "band-selective" Mott-Hubbard 
transition corresponds to the "orbital-selective" Mott 
transition in the two-orbital repulsive model [5]. But 
more than the above, what we should emphasize here 
from common features of Figs. 1 and 2 is as follows, 
(i) The metal-insulator transition between states (I) 
and (HI) is by no means "band-selective" for any ra- 
tio th/te- (ii) The position of that phase boundary on 
the plane of interactions scaled by te + th is univer- 
sal with regard to the ratio th/te. These facts indicate 
that the transition between the metallic state (I) and 
the biexciton-like insulator (HI) occurs as a result of 
the competition between the interactions and the rel- 
ative motion of electron and hole. Note that the quan- 
tity te + th is proportional to the energy of the relative 
motion. 



4. Phase diagram at arbitrary filling 

In this section, we discuss the case of arbitrary filling. 
For n ^ 1, the process for determining of the chemical 
potential /ic is added to the self-consistency cycle for 
e" and Va. Hereafter th/tc = 1 is fixed. 

Fig. 3 shows the phase diagram on the plane of 
U' and U tor n = 0.8. The Mott-Hubbard insulator 

(II) disappears immediately away from half-filling, as 
known in the single-band Hubbard model [2,3], while 
the metallic state (I) and the biexciton-like insulator 

(III) remain. In the present two-site DMFT calcula- 
tion, only the metallic state in which both Za and p" (0) 
are nonzeroes, is always obtained for U > U' within 
the present parameter region. However, it seems to ex- 
press a limitation of the two-site DMFT and the con- 
sideration in the limit of f/ — » oo actually leads the 
following results. 

In the limit of ?7 — + oo, the model (1) can be mapped 
onto a single-band attractive Hubbard model with the 
attraction —U'. From the results of the DMFT study 



3 




Fig. 3. Phase diagram in the U'-U plane for th/te = 1 and 
n = 0.8. The shaded area indicates the region speculated that 
the exciton-like insulator (local e-h pairing state) appears. 




Fig. 4. Phase diagram in the U'-n plane with th/te = 1 for the 
limiting case of large (7, referred from Ref. [7]. The interval 
area between two dashed curves is the coexistent region of the 
metallic and exciton-like insulating states. 

of this model [6,7], then, we can draw the phase dia- 
gram on the plane of U' and n in the limit of [/ — > cx), 
as shown in Fig. 4. In addition to the metallic state (I) 
and the Mott-Hubbard insulator (II) (just on n = 1 
for all values of [/'), the exciton-like insulating state 
(V) appears, in which the incoherent local e-h pairs (do 
not condense) are formed. It is worthy to note that the 
transition between the metallic (I) and exciton-like in- 
sulating (V) states is the first-order transition (except 
for n — 0, 0.5 and 1). As well as has been argued in 
Refs. [6,7], the exciton-like insulator (V) will be char- 
acterized by that Za ^ but p"(0) = (i.e., a gap 
opens at Fermi level), however which cannot be de- 
scribed within the two-site DMFT since the gap forma- 
tion {Va ~ 0) is directly connected to = through 
Eq. (4) in the present scheme. 

Let us now return to Fig 3. Based on both the result 
of Fig. 4 and the behavior of Za in the metallic state, 
it could be surmised that the exciton-like insulator ap- 
pears around the shaded region of Fig. 3, in which U' is 
comparatively large and Za of the metallic solution is 



quite small. To describe the exciton-like insulator and 
determine that phase boundary for the finite U are left 
for the future work. 

5. Discussions 

In this paper, we found the phase transitions among 
the exciton-like insulator, biexciton-like insulator, and 
the metallic state, for arbitrary filing with the use of 
the DMFT. This implies that the exciton Mott transi- 
tion can be described essentially in terms of the simple 
lattice model with only short-range interactions as well 
as the Mott-Hubbard transition. 

Finally, we discuss the relevance of our assumption 
that the e-h pairs do not condense. Although the calcu- 
lation was performed at zero temperature, we believe 
that our present results will be valid for the intermedi- 
ate temperatures, i.e., above critical temperature (Tc) 
of exciton BEC, but below temperature corresponding 
to the e-h binding energy (Eb). From simple evalua- 
tion of Tc and Eb , it can be shown that such a tem- 
perature region actually exists: consider again in the 
limit of J7 ^ 00. In the strong limit of U' , Tc can be 
estimated as of order {te + th)^ /U' [8,9]. On the other 
hand, in the low-density limit n — > 0, Eb can be esti- 
mated as of order U' [9] . Comparing these two charac- 
teristic temperatures, it can be concluded that such an 
intermediate temperature region exists even for not so 
large U' (~ te + th)- Of course, the temperature effect 
and the problem of condensation will be investigated 
by more precise calculation. 

This work is supported by CREST, JST. 
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